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A new mechanism for enhancing the hehum abundance in the blue main sequence stars of to 
Centauri and NGC 2808 is investigated. We suggest that helium enhancement was caused by the 
inhomogeneous big bang nucleosynthesis. Regions with extremely high baryon-to-photon ratios 
are assumed to be caused by the baryogenesis. Its mass scale is also assumed to be 10^ Mq. An 
example of the mechanisms to realize these two things was already proposed as the Affleck-Dine 
baryogenesis. As the baryon-to-photon ratio becomes larger, the primordial helium abundance is 
enhanced. We calculated the big bang nucleosynthesis and found that there exists a parameter 
region yielding enough helium to account for the split of the main sequence in the aforementioned 
globular clusters while keeping the abundances of other elements compatible with observations. Our 
mechanism predicts that heavy elements with the mass number of around 100 is enhanced in the 
blue main sequence stars. We estimate the timescales of diffusion of the enhanced helium and mass 
accretion in several stages after the nucleosynthesis to investigate whether these processes diminish 
the enhancement of helium. We found that the diffusion does not influence the helium content. 
A cloud with a sufficiently large baryon-to-photon ratio to account for the multiple main sequence 
collapsed immediately after the recombination. Subsequently, the cloud accreted the ambient matter 
with the normal helium content. If the star formation occurred both in the collapsed core and the 
accreted envelope, then the resultant star cluster has a double main sequence. 

PACS numbers: 98.20.Gm, 26.35.-l-c, 98.80.Ft, 13.60.Rj 



I. INTRODUCTION 

A globular cluster (GC) is a gravitationally bound sys- 
tem of stars with a spherical shape. The typical radius 
is of the order of 10 pc and its mass is about W^'^Mq. 
Stars in a GC are known to be very old and their chemi- 
cal compositions roughly represent that of the primordial 
universe. 

Among known ^ 150 GCs, uj Centauri is notorious for 
its unusual properties. Its mass is about 3 X 10^^ Mo [l| 
and it is the most massive GC in the Milky Way galaxy. 
One of the most astonishing facts about uj Centauri is 
that its dwarf stars split into two sequences the 
blue main sequence (bMS) and the red main sequence 
(rMS). 25%-35% of the dwarf stars are populated by the 
bMS [3]. Although enhanced metallicities might split a 
main sequence by lowering the surface temperature of 
metal-rich stars, it is reported that stars in the bMS are 
more metal- rich than those in the rMS in oj Centauri Q . 
Recently, it is reported that there also exists the third 
redder MS Q. This MS is suggested to be populated by 
super metal-rich stars and seems to have nothing to 
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do with their helium contents. Furthermore, observations 
by the Hubble Space Telescope discovered the triple main 
sequence in another GC, NGC 2808 j^. NGC 2808 is 
also massive, with its mass above 10^ Mq Q. It was 
found that, in either GCs, the isochrone that matches the 
properties of the bMS can be produced by assuming very 
high hehum abundances (F), 0.35 < Y < 0.45 [1 [T^- 
As other parameters do not have much effect on colors, 
it is very likely that the helium enrichment is the cause 
of the multiple main sequence. 

The origin of the enhanced helium abundances in the 
bMS stars has been a mystery. It was suggested, for 
example, that the ejecta descended from massive asymp- 
totic giant branch (AGB) stars would enhance the he- 
lium abundances of the main sequence stars and split the 
main sequence [Tl|. However, it is difficult for massive 
AGB stars to produce as much helium as F '-^ 0.35 JT^. 
What is more, the amount of helium supplied from AGB 
stars is not enough to increase the helium content of all 
stars in the bMS up to ^ 0.35 [l^- Helium enhance- 
ment by massive stars is also suggested to explain the 
MS splitting [l^. Rotating low-metal massive stars can 
enhance the helium abundance in the stellar winds but 
such stars should accelerate the helium-rich matter by 
the subsequent supernova explosions. The helium rich 
ejecta might be too energetic to be trapped in a GC and 
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could easily escape from it. These difficulties associated 
with the AGB and massive star scenarios might be re- 
lieved if the GC was originally the nucleus of a more 
massive dwarf galaxy [iS] . There were much more AGB 
stars in the dwarf galaxy and the helium supplied from 
these AGB stars might be sufficient to produce the ob- 
served number of bMS stars. As for the massive star 
scenario, a massive dwarf galaxy might be able to trap 
the helium-rich ejecta due to its deeper gravitational po- 
tential. However, there remain problems in both scenar- 
ios. It is still difficult to produce the observed fraction 
of bMS stars with the helium-rich matter or these sce- 
narios seem to require an initial mass function of stars 
extremely biased toward massive stars. In addition, no 
significant spread in metallicities of stars in NGC 2808 is 
incompatible with the observed spread in metallicities of 
stars in a dwarf galaxy. 

Most helium in the universe is believed to be pro- 
duced by the big bang nucleosynthesis (BBN). By ob- 
servations of extragalactic metal poor H II regions fl5l |. 
the mass fraction of helium produced during BBN (Yp) 
is estimated to be Yp ^ 0.25. One of the parameters 
that affect Yp is the baryon-to-photon ratio (77). Re- 
cently, results of the five-year observation of the cos- 
mic microwave background radiation (CMBR) by the 
Wilkinson Microwave Anisotropy Probe (WMAP) were 
released [l6j . The baryon density parameter is deter- 
mined to be ^bh'^ = 0.02273 ± 0.00062 where h is the 
normalized Hubble constant given by dividing the Hub- 
ble constant by 100 km s~^Mpc~^. Combining this value 
with the well-known temperature of CMBR, 2.725±0.002 
?7 is determined to be 77 ^ 6 x 10 Given this 
77, the calculation of the standard BBN yields Yp ~ 0.25, 
matching well with the observed values [1^ . In the stan- 
dard BBN theory, however, it is assumed that the BBN 
occurred homogeneously and rj has a uniform value. Al- 
though WMAP did not detect the inhomogeneity of the 
universe with scales smaller than a resolution lma.x ^ 2000 
(^'min ^ 0-1°), there remains a possibility of inhomogene- 
ity with scales under ~ 10 Mpc. Within regions smaller 
than 10 Mpc, rj might not take the standard value. Many 
calculations have been done with various values of 77 (see 
the references of ^19]). The general results are that the 
higher rj becomes, the more abundant helium will be. So, 
here comes a new candidate for the enigmatic mechanism 
of the helium enhancement in the bMS, inhomogeneous 
BBN with a high value of r;. 

We consider a baryogenesis mechanism that enhances 
the baryon density in regions with the mass scale of 
IO^Mq, a comparable amount of baryons to the GCs 
that have multiple main sequence. Though we do not 
specify this baryogenesis in more detail, the Affleck-Dine 
mechanism [20| can realize this situation right after the 
inflation era. This mechanism can enhance the baryon 
density as high as 77 ~ 1 in a region with the mass scale 
of 10® Mq [21|. The size of inhomogeneity is a free pa- 
rameter in this mechanism. The mass scale and desired 
77 indicate that the scale of necessary inhomogeneities is 



much smaller than the scales of fluctuations in the CMBR 
that can be detected by WMAP. Recent papers by Mat- 
suura et al. [l^-HJ have preformed calculations of the 
BBN with high values of 77 on the basis of the Affleck- 
Dine baryogenesis. Such an inhomogeneous BBN may be 
affected by the diffusion of neutrons. Lara [1^ showed 
that regions with enhanced 77's do not lead to abundances 
of light elements significantly deviated from those of the 
ambient medium if their size at the temperature of 100 
GK was smaller than lO'* cm. Since the size of the region 
we are concerned with in this paper is of the order of lO^'' 
cm at this temperature (see Eq. (jlOp ) and much larger 
than 10^ cm, the diffusion of neutrons cannot change the 
yields of the BBN. Still larger 77's in our scenario suppress 
the diffusion of neutrons further. In addition, this region 
extending over the event horizon even at the end of the 
BBN tolerates the diffusion of neutrons. 

In this paper, we seek a possibility that the over abun- 
dance of helium in the bMS of GCs is due to the inho- 
mogeneous BBN caused by the baryogenesis. First of 
all, we determine the parameter region that reproduces 
the helium abundance of bMS stars. The parameters 
we concerns are the baryon-to-photon ratio (77) and the 
degeneracy of electron neutrinos (^e)- £,e is defined as 

— fJ-i^^/kBT where /i^^ is the chemical potential of elec- 
tron neutrinos and is the Boltzmann constant. These 
parameters are the most effective in the sense that they 
can change the abundances of elements synthesized by 
the BBN. In Section II, we briefly summarize effects of 
the parameters 77 and during the BBN. We also sum- 
marize calculations we have done in Section II. Results 
are given in Section HI and the evolution of the helium 
enhanced region is discussed in Section IV. The conclu- 
sion is given in Section V. 

II. BIG BANG NUCLEOSYNTHESIS 

A. Free Parameters 

The BBN occurred when the temperature of the uni- 
verse was about 10^^ — 10^ K. The parameters that gov- 
ern the abundance of helium synthesized during BBN 
are 77 and for a given set of the cosmological param- 
eters. As 77 represents the density of baryons, collisions 
of baryons occur more frequently with increasing 77, re- 
sulting in larger Yp. The parameter represents lepton 
asymmetry and it can affect Yp directly through the re- 
actions 

p -I- e" < — > n + Vo, (1) 

which determines the number of neutrons and protons 
before BBN. At the equilibrium of these reactions, the 
number ratio of neutrons to protons is expressed by a 
function of the temperature T and as 

^ = exp(-^-ee), (2) 
Up \ k^T J 
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Table. I: Nuclei included in the small reaction network. Z is 
the name of nuclei and A is the mass number of Z included 
in the network, 'n' refers to neutron. 



where Arac^ is the difference of the rest energies between 
a proton and a neutron. As most of neutrons that had 
existed before the BBN era were brought into hehum nu- 
clei, a negative enhances helium synthesis. We treat 
as a free parameter restricted from observations [1^ and 
assume /3-equilibrium at the beginning of calculations. 
Allowed range with 2a CL is 

-0.01 < < 0.1. (3) 

Note that the Affleck-Dine mechanism prefers negative 
with the absolute value of the order of t] because t] is 
supposed to be positive W\ . 



B. Numerical Calculations 

Because the horizon scale right after the BBN is much 
smaller than lO^M© in mass, the effect of diffusion is 
negligible. Thus, we calculate the BBN with high r\ in 
the same way as the homogeneous BBN ignoring inho- 
mogeneities inside the region and adopt the same method 
of numerical calculations as that of the ordinary homo- 
geneous BBN. We first calculated BBN with the reaction 
network including 61 nuclei shown in Table. ID The nu- 
cleosynthesis with a larger reaction network with nuclei 
in Table. |ll]is also calculated to see how much heavy el- 
ements are produced. Heavy element production during 
BBN with high r\ is discussed in [23 - l24| . To calculate 
BBN, we used the Kawano code introduced in [2^, with 
revised reaction rates taking from [29] whenever they are 
available. Reaction rates of heavy elements are taken 
from [30| . 



III. RESULTS 

We calculate the BBN with various parameters. The 
calculations are started from the time when the tempera- 
ture is 10^^ K and we obtained the abundances at lO*" K. 
The nuclear statistical equilibrium and the /^-equilibrium 
are assumed at the initial state. 
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Table. H: Nuclei included in the large reaction network. See 
the caption of Table. Ufor details. 

A. Lower Limits 

To locate a star on the bMS in w Ccntauri by enhancing 
the helium content, the helium fraction should be as high 
as 0.35 Q. Therefore we require 77 to be greater than 77b 
that reproduces = 0.35. To estimate the amount of 
helium produced during BBN, only the reactions of light 
elements are essential. For this reason, we use the small 
reaction network to estimate Yp. 

The results of calculations are shown in Fig. [TJ We 
performed the calculation with three different fe! = 0, 
0.01, and —0.1. The horizontal line indicates = 0.35. 
For = 0, 77 yielding Yp > 0.35 is in the region 

77 > 3 X 10"^ (4) 

Thus 77b = 3 X 10"^ For = -0.01 and = 0.1, 
77b = 2 X 10^^ and 77b = 2 x 10^^, respectively. 77 above 
these values are required to explain the bMS stars. 

B. Upper Limits 

As previously predicted by Matsuura et al. (22| - [23 |. 
heavy elements are produced during BBN with high 77. 
The amount of heavy elements produced during BBN 
gives the upper limit of 77. Too large 77 might overproduce 
heavy elements incompatible with the metallicity of bMS, 
which is at least a factor of a few smaller than the solar 
abundances. As nuclei with the mass number [ A) of ~ 80 
or greater are important for this restriction, we use a 
large reaction network for this purpose. For calculations 
of the large reaction network, we only consider the case 
of = 0. 

In Fig. [21 we show the results of BBN calculations for 
some values of 77. Between 77 = 1 x 10"^ and 77 = 2 x 10~^, 
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Fig. 1: Results of Yp for some sets of ri and ^e- The horizontal 
line at the middle represents Yp = 0.35. 



Fig. 2: The number fractions of the produced elements rel- 
ative to the solar abundances. All the calculations are done 
with ~ 0. 



there is a considerable jump in the amount of heavy el- 
ements with A ^ 100. Fig. 7 of Matsuura et al. ^ 
shows that rj = 10^* yields a significant amount of heavy 
elements with A ~ 150, comparable to the solar abun- 
dances. With higher rj, the abundances of heavy ele- 
ments continue to increase. Observations show that the 
amounts of heavy elements are at least a factor of a few 
less than the solar abundances in the relevant GC _3l|. 
Therefore the upper limit of 77 is somewhere between 
3 X 10-5 and 10"^. 

In Fig. [21 some elements A ^ 100 are very close to the 
expected abundances of the GC stars. Thus, it is a good 
touchstone of our model to check the abundances of the 
elements A ^ 100 in the dwarfs in bMS. 



A. Diffusion 

Helium ions in the region with a higher l^p might dif- 
fuse outward the surrounding region with lower helium 
contents. This effect takes place after the horizon scale 
exceeds the size of the region, that is, after the tem- 
perature becomes lower than 10* — 10^ K. The diffusion 
equation for helium is written as 



dn 



He 



dt 



(5) 



where tt-hc is the number density of helium and D is the 
diffusion coefficient of helium with electrons and its mag- 
nitude is estimated as 



D = 



(6) 



IV. EVOLUTION OF HELIUM ENHANCED 
REGION 



Even if a primordial gas with high Y-p was once formed 
at the time of BBN, it might be mixed with other gases 
in the outer regions and end up with a lower fraction 
of helium. We briefly show that the effect of such mix- 
ing would not affect the region with high helium abun- 
dance. Here, we do not only consider effects of diffusion 
and accretion, but also take into account self-gravity as 
a competing process. 



where A is the mean free path of helium ions and v = 
\J ik^T / Am]i is the thermal speed of helium ions. Here 
the constant toh denotes the mass of hydrogen. The time 
scale tdiff of the diffusion is approximated as 



idiff — 



d In riHo 



dt 



D 

where R is the size of the region. 
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Before Recombination 



2. After Recombination 



Before the recombination of electrons, photons and 
matters had been coupled and had the same tempera- 
ture. The number density of helium during this epoch is 
proportional to the cube of the temperature and written 
as 



To 



(8) 



Here To = 2.73 K is the present temperature of CMBR 
and n^o is the present number density of photons and 
expressed as n^o = 2C{3)/tt'^ {k-BTo/Hc)^ — 416 cm^'"^ 
where h is the reduced Planck constant. For the case 
of 7/ = 10-* and Fp = 0.4, 



nHc = 2 X lO-^T^ cm 



(9) 



Thus, the size R is expressed as a function of the tem- 
perature 



R 



\ 167^mH^^He 



1/3 



5 X 10^1 
f 



(10) 



Next, we estimate D before the recombination. As- 
suming that the universe has no net charge, the number 
density of electrons is written as 



Tie = + 2n 



He 



= 1- 



T 

2 J 

2 X lO^^T^ cm"3 



The last value is, again, for the case of rj 

8 X lO^T^ 



Yp = 0.4. From [3; 



A 



4nee4 In A 



Uf, In A 



cm. 



(11) 
(12) 

(13) 
10-'' and 

(14) 



where 



A = 



4e2 



1/2 



= 7 X 10' 



1/2 



Thus, 



D = 



Xv 



X 10 



11 



cm^ sec ^. 



(15) 



(16) 



From R and D derived above, t^is is estimated to be 
3 X 10^1 



tdiS 



Vio^K ) 



sec - 9 X 10^° ( 1 sec. (17) 



This diffusion time scale is so large that the effect of 
diffusion is negligible. 



Matters are assumed to completely decouple from ra- 
diation after the recombination. Thus, njjc and R after 
the recombination are denoted as functions of the gas 
temperature T, 



"-Ho = — 7?%0 ^ 



3/2 



30^,3/2 



and 



„ 7 X 1019 



(18) 



(19) 



Here the temperature Tro when free electrons recombine 
with hydrogen ions is about 5,500 K due to the higher 
density. We have assumed that all of helium and hy- 
drogen atoms are neutral in this stage. Therefore the 
cross section for helium atoms becomes of the order of 
^ 10~i^ cm^. Accordingly, D becomes 



D = 



3 X 10*^ 



T 



cm^sec ^. 



(20) 



The diffusion timescale is obtained as a constant value. 



(21) 



tdiff = 2 X 10^^ sec. 
This diffusion time is also too large. 

B. Accretion 



Since the region with enhanced helium content has a 
much higher density than the surrounding region, the 
region accretes matters from the outside, which might 
dilute the helium fraction. For simplicity, the timescale 
of this process is estimated assuming the Bondi accretion 
(33} . The time scale of the Bondi accretion is inferred 
from the relation 



= 4.g^m4^A(7), 



M 



!(oo) 



(22) 



where M is the total mass of the region (lO^M© in this 
case), p(oo) and Cs(oo) are the density and the sound 
speed at infinity. The eigen value of the steady state so- 
lution A(7) is a function of the adiabatic index 7 written 
as 



A(7) 



(5-37) 



=5/3 



(7^1) 
(7 = f )• 



(23) 



The values of p{oo) and Cs(oo) are taken from those in 
the surrounding region with Yp = 0.25 and 77 = 6 x 10~^'^. 



6 



They change as the universe expands and can be written 
as functions of temperature, 



p(oo) 



(24) 



2 X 10"32r3 gcm- 



(before recombination) 
(25) 

3 X 10 27y3/2 g(.j^ 3 (af|;gj- recombination), 

where T after the recombination denotes the gas temper- 
ature. Similarly, Cs(oo) is expressed as 



Cs(oo) 



3 1- 



2x10^" 



3x 103 



cm sec ^ (before) 



(26) 



(27) 



10''\/r cmsec-i (after), 



where fi is the mean molecular weight and ub is the ra- 
diation energy constant cb = Sir^kg/lbh^C^ . Thus, the 
accretion time scale iacc becomes 



3 X 10^6 (l + 3^^^^^ T-3 sec (before) 



_ M _ 

tacc — ~T~ 

M 



the other hand, the Jeans mass in this region reduces to 
2 X 10"^ Mq immediately after the recombination. The he- 
lium rich region collapses on the time scale of^ 1.5x10^^ 
sec, much shorter than the accretion time scale at the 
temperature of T^c — 5,500 K, i.e., 3.5 x 10^^ sec (see 
equation (1281) for before recombination). This preceding 
collapse might avoid overall mixing of the helium rich 
matter with the accreted ambient matter. The resultant 
cloud consists of helium rich core surrounded by a pri- 
mordial gas with the normal helium content. 

Though there is no widely accepted scenario for the 
formation of a GC, we raise possible routes to realize 
GCs with multiple main sequence. Subsequent star for- 
mation both in the core and envelope could lead to a 
star cluster with double main sequence resembling uj Cen- 
tauri. To reproduce a GC like NGC 2808, two such clouds 
with different Yp's in helium rich cores are needed to col- 
lide and trigger star formations in each of regions with 
three different Yp's. The collision with a high velocity 
removes stars that have formed before. This may also re- 
duce the star-to-star scatter in metallicities. As a result, 
stars in the thus formed cluster have main sequences with 
three different helium contents and have similar metal- 
licities. These features are observed for stars in NGC 
2808. Cloud-cloud collisions have been raised as a pos- 
sible mechanism for the globular cluster formation [34| . 
Young stellar clusters found in merging gala xies [35| and 
near the center of the Milky Way galaxy [3^] support this 
mechanism. 



(28) 



6 X 10^^ sec 



(after). 



This result means that the time scale of accretion be- 
fore the recombination is long enough to avoid dilution 
of helium, while the time scale after the recombination 
becomes significantly shorter than the Hubble time at 
that epoch. This indicates that the average helium con- 
tent reduces considerably. However, if we take the effect 
of self-gravity into account, star formation might precede 
the dilution of helium depending on the effectiveness of 
turbulent mixing. 



C. Star Formation 

The Jeans length Aj of the system before the recom- 
bination is 



A,, = 



and the Jeans mass Mj is 
4 



1027^-3/2 



1 + ^ 



(29) 



- x3 8x1021Mq 
3 {T + 3xmf^ 



(30) 



Therefore as long as the matter is ionized the mass of 
the helium rich region never exceeds the Jeans mass. On 



V. CONCLUSIONS 

We calculated the BBN for several values of rj and 
and looked for the parameter region where the helium en- 
richment in the GCs is reproduced. We set lower limits 
of to be 77b = 3 X 10"^ for = 0, 77b = 2 x 10^^ 
for $e = -0.01, and 77b = 2 x 10"^ for = 0.1 
. The upper limits come from metal abundances and 
the allowed parameter region for the case of = is 
3 X 10~^ < 77 < 10~^. If an inhoniogeneity caused by the 
Affleck-Dine baryogenesis made a region with 77 and in 
these criterions, the helium abundance is enhanced and 
the region can end up with GCs, showing multiple main 
sequences. Our model predicts that the abundances of 
heavy elements with A ^ 100 is enhanced in the bMS 
stars. Stars in our model are expected to have more con- 
centrated spatial distribution of the bMS stars than the 
rMS as was recently reported for oj Centauri because 
the bMS stars are formed in the central part of the he- 
lium enhanced region while the rMS stars are formed in 
the outer envelope in our model. 

A problem of the Affleck-Dine baryogenesis is that it 
can be adopted to antimatters as well as ordinary mat- 
ters and can create regions with enhanced antimatters. 
As there are no observations of such regions with an- 
timatters, there should be some mechanism that stops 
the generation of antimatters or there must be another 
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baryogenesis mechanism to exclusively enhance ordinary 
matters up to ~ 3 x 10~^ in small regions. 
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